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We study the dynamics of the Mott insulator-superfluid 
quantum phase transition in a periodic ID array of Joseph- 
son junctions. We show that crossing the critical point at a 
finite rate with a quench time tq induces finite quantum fluc- 
tuations of the current around the loop proportional to Tq^^^ . 
This scaling could be experimentally verified with in array of 
weakly coupled Bose-Einstein condensates or superconduct- 
ing grains. 
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Thermal fluctuations critically slow down close to a 
continuous symmetry breaking thermodynamic phase 
transition. Both the relaxation time t and the corre- 
lation length ^ diverge when a dimensionless parameter 
e, which measures the distance from the critical point, 
tends to zero: ^ ~ le]""^ and r ~ kl"^", with z and v 
critical exponents. In the infinitesimally slow transition 
the symmetry broken phase is entered in a state of equi- 
librium with a fully ordered complex order parameter. 
However, in a quench occuring at a finite rate transition 
the critical slowing down (r — > cx) when e — > 0) implies 
that the system goes out of equilibrium some time before 
the transition. As a result the complex order parame- 
ter which emerges after the transition will assume differ- 
ent random phases in different domains of space. This 
leads to formation of topological defects such as vortices 
wherever the circulation of the phase around a closed 
loop happens to be nonzero. This process is known as 
Kibble-Zurek mechanism (KZM) ||-||. The size of the 
domains follows from a simple argument Q : Close to the 
transition one can linearize the time dependence of e as 
e{t) « t/TQ. The transition rate is r{t) = e/e = l/\t\. 
The system goes out of equilibrium at the time —t be- 
fore the transition when r(—t) = T~^{~t). After that 
time the state of the system essentially does not change 
until +t, when the rate r{t) becomes again equal to the 



relaxation rate, r{-\-i) = t At i = 



fluc- 



tuations of the order parameter with wavelengths longer 
than i = ^{i) = Tq^<^^+'''^ begin to grow exponentially 
while short wavelength fluctuations remain unchanged. 
The size of the correlated domains is given by ^ and the 
density of vortices by 1 /^^ . 

Theoretical and experimental studies of KZM have 



concentrated so far on thermal continuous phase transi- 
tions with the dynamics of the order parameter governed 
by an effective irreversible time-dependent Ginzurg- 
Landau theory. Most attention has been devoted to the 
normal-superfluid transition in '^He superconductors 
|5|, and, more recently, dilute Bose-Einstein condensates 
|7|. To date, truly microscopic quantum approaches have 
been too complicated to extract useful predictions. 

A quench-induced quantum phase transition (QPT) at 
temperature T = must be treated in a microscopic way. 
It is a common wisdom that some properties of a quan- 
tum transition can be obtained by an exact map from 
a thermodynamic transition Q|: the correlation length 
in the ground state of the quantum system scales like 
^ ^ \e\~'^ and the gap between the ground state and the 
first excited state like dE ^ \e\^^'^ . To study the dy- 
namics of a quantum transition, which drives the system 
out of its ground state, we also need information about 
correlations in excited states, which is not provided by 
this map. As we will see below, the essence of the KZM, 
which is the competition between the transition rate and 
the timescale on which the system can react, is applicable 
to quantum transitions. However, the quantum scenario 
and in particular interpretation of its results differ from 
the thermal case. The main reason is the reversibility 
of the quantum dynamics, as opposed to the dynamical 
irreversibility of the thermal critical dynamics. 

In this Letter, we study the appearance of a 
nonzero current while the system is undergoing a (zero- 
temperature) quantum phase transition. We develop a 
microscopic, dynamical theory and suggest experiments 
to test our predictions. 

Josephson junction arrays. A prototype system 
displaying a continuous QPT is an array of mesoscopic 
Josephson junctions (JJ). Superconducting J J networks 
are now available thanks to important advances in tech- 
niques for ultra-small superconducting grains although 
existence of gauge fields may complicate the analysis [|j . 
On the other hand, manipulation of optically trapped 
Bose-Einstein condensates promises the observation of 
QPT in neutral quantum fluids [^0|jl^]. In such systems, 
a superconducting/superfluid-insulator phase transition 
is driven by the competition between two physical magni- 
tudes: the Josephson coupling energy Ej, which governs 
the tunneling through the intra- well barriers, and the on- 
site interparticle interaction energy Ec- When they are 
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comparable, there is a competition between long-range 
order (which is favoured by Ej) and localization (induced 
by Ec)- When the latter prevails, no net current can flow 
through the junctions. 

An effective Hamiltonian describing the quantum dy- 
namics of an array of J J's is given by the Quantum Phase 
Model (QPM) 

;=i {k,i) 

with the last sum running over the nearest-neighbour 
sites. For our purposes, it is enough to consider the 
case of a diagonal charging energy matrix E^ indepen- 
dent of the site indices. The phase 0; and the num- 
ber of atoms hi in each site of the junction are (with 
some important caveats fl^ ) non-commuting conjugate 
observables = i, and, in the (/)-representation, 

hi = n ~ ^^7! 4'i = 4*1 ■ Therefore, (jji and n; play the 
role of coordinate and conjugate momentum, and sat- 
isfy the Heisenberg uncertainty relation The insulator 
phase is characterized by large quantum phase fluctu- 
ations in the ground state, which destroy the long-range 
order among sites. The ID Hamiltonian Eq.(^ exhibits 
a continuous Mott phase transition at the critical value 
Gc = Ej/Ec = 0.617 0. In what follows we study the 
microscopic dynamics of quantum phase transitions in 
a one-dimensional chain with periodic boundary condi- 
tions. 

The dynamics is governed by the QPM Hamiltonian 
Eq.(n]), and satisfies a Schrodinger equation 




with G ~ Ej/Ec and the time rescaled a.s t ^ Ect. This 
equation is valid for a large number of atoms per site, 
n > 1, and for G <^ ri^ 0. The ground state be- 
low the transition [G < Gc) is close to the Fock state 
\n, n,n, . . .) with all sites occupied by the same number 
of atoms n. In the phase representation this state is de- 
scribed by a uniform wavefunction, ^'(0/) = const, where 
all phase differences between nearest neighbor sites have 
the maximal dispersion of Acjj ^ 1. Above the transition 
(G > Gc) the ground state is a number squeezed state 
which continuously tends to a Fock state when G — > G+ , 
and to a coherent state for G ^ n^. Therefore, when 
G ^ Gc one can describe the low energy part of the 
spectrum of Eq.(^ in a harmonic approximation |^5| 

^ 1=1 ^ (kd) 

This equation is diagonalized by normal modes num- 
bered by momentum n £ {—Ns + 1, . . . , -|-A^s}, 5* = 



n^^A'(^A')- There is one zero mode $0 ^ ^i4>i- 
All other modes $^ have nonzero frequencies ^7^G[« 
Qi/2(27^^ for small ^], which scale as G^^^, 

To implement quench with a quench timescale tq we 
linearly ramp the control parameter G in Eq.(^ as 

G{t) = - . (5) 

The time t runs from to rgGmax with Gmax S> Gc. 

Adiabatic transition. In the limit tq ^ 00 the tran- 
sition is adiabatic and the system state adiabatically fol- 
lows its ground state from the Fock state |n, . . . ,rt) at 
G = to a coherent state with all the nearest neighbor 
phase differences close to zero, A<^ « 0, for Gmax ^ Gc- 
The ground state at Gmax ^ Gc is well described by the 
ground state of the harmonic oscillators (Q). The disper- 
sion of phases in any harmonic oscillator ground state at 
Gmax is proportional to Gmax^- The final dispersion of 
phase differences in the adiabatic transition is 

Mrc^^oo ^ G^J^' . (6) 

Instantaneous transition. In the opposite limit of 
an istantaneous transition, when tq — s- 0, the system 
is not able to adjust its quantum state to the changing 
hopping rate G(i), and it remains in the initial Fock state 
till Gmax- The final dispersion of the phase differences in 
the uniform Fock state, '^{(t>i) = const, is 

A^rQ^O ^ 1 . (7) 

For a generic quench the evolution of the system is 
approximately adiabatic when G is far from Gc, and im- 
pulse when G is close to Gc', compare Fig.l for a 3-site 
periodic lattice. This is similar to the thermal case 
As G increases starting from the system initially follows 
its ground state, which is an incoherent superposition of 
Fock states, until a G_ < Gc when the transition ceases 
to be adiabatic. Between G_ and a certain G > Gc the 
evolution is impulse: to a first approximation the state 
of the system does not change. The system arrives at 
G in an incoherent superposition of Fock states with a 
dispersion ~ 1 as in the initial Fock state. After G 
the transition again becomes adiabatic. 

Fast transitions with {Ec)tq ^ 1. Let us esti- 
mate G for a dimensionless tq <C 1 when we expect that 
G ^ Gc and we can use the harmonic approximation 
(^,^ for G ~ G. For the linear quench (||) the peri- 
ods of the oscillators decay like G{t)^^^^ = ^rg/iGmax- 
The transition timescale is G/G = t. This rate becomes 
comparable to the periods of the oscillators at the time t 
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when G ^1"^ ~ G/G: compare Fig.l for a 3-site periodic 
lattice. A solution of this equation gives 



rQ«l 



r]l\ or i?-/3,V3) ^ 



Gtq-^I = G(<rQ<Sl) — '''q ■ 

As anticipated, when tq <C 1 we have Gt-q^^i 3> Gc 

/ 



(8) 
(9) 




FIG. 1. The inverse of the gap T = as a function 

of G for a 3-site periodic Hubbard model with n = & atoms 
per site (dots). The solid line is the best G~^^'^ fit. The 
dashed lines show the transition time G/G = t — Grg/Gmax 
for Gmax = 10 and tq = 4.5, 2.1, 0.9 (from left to right). The 
crossings between 1/dE and the dashed lines define G for the 
different tq. 

We can use the fact that GrQ<i to estimate the fi- 
nal dispersion A</) at Gmax- At GrQ<i the dispersion 
is A(/) ~ 1. For G > GrQ<i the evolution of the har- 
monic oscillators is adiabatic. There is no mixing be- 
tween the eigenstates of any oscillator /i (Q). The dis- 
persion of the phase in any given eigenstate of the 
oscillator (^ scales as G^^/"'. In the adiabatic evolu- 
tion for G(t) > GrQ<gi the dispersion of <i>^ shrinks like 
[GrQ<gi/G(t)]^/^. shrinks in the same way. The final 
dispersion at Gmax is ~ [GrQ«i/G„iax]^/'' or 



A</), 



TQ<1 



,^/«G-V^( or E:''\^"'G-Ji') . (10) 



In a ID periodic lattice this l^(j)TQ<^i translates into a 
dispersion of the angular momentum L — z^j(a|a/+i — 
Oq+iO-i) through the formula 



AL 



'Ns n Act) 



(11) 



This Ai may translate into dispersion of the winding 
number \/NsA(j) when the atoms after the quench are 
forced to condense. Eq.([l0|) is valid when the predicted 

" 2/3 . 

Gtq<i Tq is much less than the final Gmax- This 

condition is satisfied when tq 3> Gmax . In fact for the 
slowest /X = ±1 modes to become adiabatic before the 



maximal possible G,, 



we need tq ^ n/Ns 



Figures 2 and 3 show results of numerical simulations 
of the 3-site periodic lattice with a total of 300 atoms. 



In Fig. 2 we show that A(j> remains the same as in the 
initial Fock state for G < GrQ<i. In Figure 3 we verify 

1/3 

the scaling t ~ Tq . We also performed similar simu- 
lations for a total of 120 (commensurate) and 120+1 
atoms (non-commensurate density of atoms). Plots like 
those in Fig. 2, 3 for 120 and 121 atoms arc impossible to 
distinguish which shows that an extra noncommensurate 
atom does not make any difference for fast transitions. 




FIG. 2. The normalized dispersion dl — AL/ ALtq^o 
as a function of S = G(t)/Gniax. We show results of ex- 
act numerical simulations of the 3-site Hubbard model for 7 
different tq = 0.016, 0.008, 0.004, 0.002, 0.001, 0.0005, 0.00025 
(from left to right) with Gmax ~ 1000 and n = 100 atoms per 
site, dl deviates from 1 (and hits the bottom of the figure) at 



G. 



-2/3. 



compare Fig. 2. 




FIG. 3. The normalized dispersion dl — AL/ ALtq~*o as a 



function of s = G(t)/G with G ~ t^ 



-2/3 



This figure displays 



all the seven plots from Fig. 2. The plots in Fig. 2 are in the 
top left corner of this figure. The fact that the seven plots sit 
on top of each other proves that the evolution of AL depends 
on time t through the combination t/t with t t^ 
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Slow transitions with {Ec)tq ^ 1. So far we con- 
centrated on TQ <C 1 with G well in the harmonic regime. 
In slower transitions G is close to Gc where the harmonic 
approximation cannot be applied. The slow transitions 
(tq <C 1) are probing the critical behaviour close to Ga- 
in the perioric 3-site array there is no phase transition 
but, rather, a crossover: the gap dE is minimal at the 
crossover point Gc « 0.5 but it does not vanish there. 
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compare Fig.l. For tq ^ 1 the transition is adiabatic 
and the system follows its ground state. The phase dif- 
ference dispersion at Gmax ^ Gc is the dispersion in the 
ground state of the harmonic oscillators (^). A thermal 
crossover transition was studied in |19|. 

In a large 2 D array with a commensurate density of 
atoms there is a phase transition and not a crossover 
|@,|0|. The gap dE vanishes at G = Gc- Close to 
the critical point, G ~ Gc, the energy gap scales as 
dE ~ (G — GcY , where v = 2/3 is an exact renormaliza- 
tion group critical exponent. Note that in the harmonic 
regime for G 3> Gc we have dE ^ G^l'^ which is consis- 
tent with a mean- field value oiv = 1/2. With a linearized 
G{t) -Gc = ^ the gap Ai5 ~ (G - GJ^/s = (t/rQf/^ 

becomes equal to the quench rate r{t) — G/{G — Gc) = 
1/t at G ~ Gc -f . Before G the system is in an in- 

coherent superposition of Fock states. After G the evolu- 
tion is adiabatic, and with increasing G the phase disper- 
sion shrinks together with the phase widths of the system 
eigenstates. The dispersion measured at Gmax ^ Gc, 
where we can use the harmonic approximation, is 



(Gc + 



(Be 



1/4 



rQ>l 



G 



1/4 



G 



1/4 



(12) 



This formula is consistent with Eq.(^ because Gc = 
0(1). Here we again use the fact that the phase disper- 
sion of harmonic oscillator eigenstates shrinks like G^^/^. 
The critical behaviour is realized in the {Ec^^)tq^ term. 

Concluding remarks. The state of the system after 
the transition does not have definite angular momentum 
or definite phase differences between lattice sites. That 
would be the case in the thermal KZM, where Acj) would 
describe the dispersion in an ensemble of different pos- 
sible classical outcomes. In a quantum phase transition 
at zero temperature the state of the system after the 
transition is in a coherent superposition of states with 
different A(j>. Either a measurement or decoherence ||2^] 
are needed to convert this coherent superposition into a 
mixtures of states, each with definite current. 

Another difference with respect to the thermal KZM 
is due to the reversibility of the quantum dynamics. In 
the thermal case the characteristic lenghtscale ^ is frozen 
after the symmetry breaking transition is completed. ^ is 
a permanent record of the transition rate tq: manipula- 
tions with e do not change the winding number as long as 
the system remains in the symmetry broken phase with 
e > 0. In contrast, even after the diabatic quantum tran- 
sition is completed one can change A(/)^q . Adiabatic vari- 
ations of G away from Gmax a-re accompanied by changes 
in A0, A(/) = A0,^(Gmax/G)i/4. 

In conclusion, we have predicted the phase dispersion 
after a diabatic insulator-superfluid quantum phase tran- 
sition in an array of Josephson junctions. This theory is 
a quantum counterpart of the Kibble-Zurek mechanism 



for topological defect formation in classical thermal phase 
transitions. Our predictions can be tested experimentally 
in ID or 2D superconducting J J arrays. The possibility of 
tuning the Josephson coupling energy in a superconduct- 
ing J J has been demonstrated recently iQ. ID (or 3D 
p2[) J J arrays realized with Bose-Einstein condensates 
can be, perhaps even more easily, accurately tailored 
|p^,pI|jT7|]l|]. For the case of a ID array our Eq.(|l^) pre- 
dicts a dispersion of the phases A^ ~ Gmax^Tg™ that 
can be directly measured by the interference techniques 
which have been used in the experiments ||l^,|ll|,0 . 
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the DOE and the Cofinanziamento MURST. 
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